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ABSTRACT 
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The  Harry  Dym  equation,  which  is  related  to  the  classical  string  problem, 
is  derived  in  three  different  ways.  An  implicit  cusp  solitary  wave  solution 
is  constructed  via  a  simple  direct  method.  The  existing  connections  between 
the  Harry  Dym  and  the  Korteweg-de  Vries  equations  are  uniformised  and 
simplified,  and  transformations  between  their  respective  solutions  are  carried 
out  explicitly.  Whenever  possible,  physical  insights  are  provided.  /  u  j 
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1.  Introduction  mL 

In  trying  to  discover  possible  nonlinear  PDEs  that  could  be  solved  using  our  real 
exponential  approach  (Hereman  et  al  1985, 1986),  we  stumbled  upon  the  Harry  Dym  (HD) 
equation  rt  =  r^rsz  that  has  nonlinearity  and  dispersion  coupled  together.  This  prototype 
of  an  evolution  equation,  that  admits  a  cusp  solitary  wave  solution  (Wadati  et  al  1980), 
appears  in  many  disguises,  namely,  r*  =  (1  —  r)®rsx,  u  =  (r“a)si  and  (r^)t  =  (r“^)32. 
The  first  one  is  occasionally  referred  to  as  the  cusp-soliton  equation  (Kawamoto  1984a,b, 
1985a).  True,  it  failed  to  fit  in  our  scheme  but  it  certainly  kindled  our  interest.  According  to 
some  of  the  early  references  (Kruskal  1975,  Magri  1978,  Wadati  et  al  1979,  1980,  Sabatier 
1979a,b,  Dijkhuis  and  Drohm  1979,  Case  1982,  Yi-Shen  1982,  Calogero  and  Degasperis 
1982)  the  origin  of  the  HD  lies  in  various  private  communications  with  Harry  Dym.  In 
search  of  truth,  we  contacted  Dr.  Harry  Dym  who  replied,  the  spring  of  1974,  ... 
Martin  Kruskal  delivered  a  few  lectures  on  the  isospectral  theory.  ...  Motivated  by  these 
lectures,  I  developed  some  analogues  for  the  string  equation.  The  HD  equation,  as  Martin 
later  termed  it,  was  one  of  the  outcomes.”  Their  collaboration  resulted  in  developing  a 
draft  of  "a  fairly  complete  theory”  whidi  is  ”gathering  dust”  in  Dr.  Kruskal’s  office. 

fkom  the  fact  files,  the  HD  is  a  completely  integrable  nonlinear  evolution  equation 
(Wadati  et  al  1979, 1980)  which  can  be  solved  by  the  Inverse  Scattering  Transform  (1ST). 
Dijkhuis  and  Drohm  (1979)  and  Calogero  and  Degasperis  (1982)  discuss  the  HD  equation 
as  a  special  case  of  a  new  broad  class  of  nonlinear  PDEs  tractable  by  1ST.  The  HD  has  a  bi- 
Hamiltonian  structure  (Magri  1978,  Case  1982,  Olver  1986);  it  possesses  an  infinite  number 
of  conservation  laws  and  infinitely  many  symmetries  (Magri  1978,  Ibragimov  1985,  Olver 
1986);  and  reciprocal  Backlund  transformations  (Rogers  and  Nucci  1986,  Nucci  1988). 
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However  the  HD  equation  does  not  possess  the  Punlev4  property  (Weiss  1983, 1986,  Steeb 
and  Louw  1987,  Fokas  1987,  Hereman  and  Van  den  Bulck  1988)  indicating  that  the  Painleve 
property  is  at  most  sufficient,  but  not  neccesary,  for  integrability. 

Yi-Shen  (1982)  and  Sabatier  (1979a,  b,  1980)  recall  the  derivation  of  the  HD  equation. 
The  first  author  starts  from  the  isospectral  equation  modelling  the  classical  string  problem 
with  a  varying  elastic  constant.  Sabatier  independently  rediscovers  the  HD  and  its  general¬ 
isations  within  the  Lax  formalism  (Lax  1968).  In  order  to  make  this  paper  a  self-contained 
tutorial  of  the  HD  equation,  we  shall  first  review  and  simplify  these  calculations,  and  also 
include  a  heuristic  derivation  of  the  HD  to  provide  some  physical  motivation. 

Amongst  the  first  researchers  to  study  the  connections  between  the  HD  and  the  KdV 
equations  we  mention  Ibragimov  (1981,  1985)  Calogero  and  Degasperis  (1982)  and  Weiss 
(1983, 1986).  Later  on  Kawamoto  (1985b),  Rogers  and  Nucci  (1986)  and  Steeb  and  Louw 
(1987)  investigated  the  links  to  the  mKdV  equation.  Underlying  all  these  connections, 
there  exists  a  rigorous  though  involved  Lie-Backlund  algebra.  Excellent  references  on  this 
approach  are  in  books  by  Ibragimov  (1985)  and  Olver  (1986).  In  keeping  with  a  tutorial 
point  of  view,  we  therefore  recognise  the  need  for  uniformisation  and  simplification  of  these 
connections  in  order  to  provide  a  better  insight.  Abo,  as  part  of  our  simplifying  efforts,  we 
rederive  the  implicit  solution  to  the  HD  equation  in  the  simplest  possible  way,  without  the 
aid  of  heavy  machinery  like  the  1ST.  We  obtain  a  particular  exact  solution  which  cannot, 
however,  be  expressed  in  closed  form  owing  to  the  presence  of  a  transcendental  phase. 
Although  the  links  between  the  equations  have  been  extensively  studied,  limited  attention 
has  been  psud  to  connections  between  the  various  solutions.  In  our  opinion  this  is  vital, 
since  a  customary  technique  in  nonlinear  science  b  to  generate  solutions  to  a  nonlinear 
PDE,  from  known  solutions  of  another. 

The  organbation  of  thb  paper  b  as  follows.  In  Section  2,  we  first  formally  derive 
the  HD  from  the  string  problem  with  a  varying  elastic  constant  and  where  the  eigenvalue 
A  b  constant  w.r.t.  a  parameter  t.  We  next  employ  the  Lax  operator  technique  for  this 
isospectral  eigenvalue  problem.  Lastly,  we  retrieve  the  HD  by  a  heuristic  method.  This 
involves  the  derivation  of  the  linear  part  of  the  PDE  from  the  known  dbpersion  relation, 
and  suitable  modification  of  the  coefficients  to  accoimt  for  nonlinear  effects. 

As  stated  earlier,  the  HD  b  conventionally  solved  using  1ST.  We  present,  in  Section 
3,  a  novel  direct  integration  method  to  construct  the  implicit  cusp-type  single  solitary 
wave  solution  of  the  HD  equation.  Guided  by  an  a  pnort  knowledge  of  the  final  result,  we 
assume  a  form  of  the  solution  which  b  inherently  implicit.  Mathematically  speaking,  this 
b  achieved  by  a  change  of  variabb  that  depends  on  the  solution  of  the  equation  itself.  As 
will  be  explained  later,  the  introduction  of  thb  type  of  variable  b  commensurate  with  the 
physical  basb  for  implicit  solutions  of  other  nonlinear  PDEs,  vb.,  the  nonlinear  kinematic 
equation  (for  other  examples  see  e.g.  Wadati  et  al  1979,  1980  and  Kawamoto  1984a,  b. 
1985a). 

Now,  the  HD  equation  has  intriguing  links  with  other  nonlinear  evolution  equations, 
vb.,  the  ubiquitous  KdV,  and  a  Liouville-type  equation  (Hereman  and  Banerjee  1988). 


In  Section  4,  we  therefore  summarise  some  of  the  existing  links  between  the  HD  and 
the  KdV  by  transforming  one  into  the  other  and  by  providing  the  connections  between 
the  corresponding  eigenvalue  problems.  Scattered  work  done  by  Ibragimov  (1981,  1985), 
Weiss  (1983,  1986)  and  Kawamoto  (1985b)  have  contributed  to  the  overall  picture  we 
draw  in  this  section.  The  procedures  we  employ  comprise  the  Backlund  transformation 
method,  a  technique  involving  the  Schvoarzian  derivative,  and  the  Cole-Hopf  and  Miura 
transformations.  Every  procedure  above  needs  to  be  augmented  by  either  an  explicit-to- 
implicit  or  an  implicit-to-explkit  transformation.  It  is  thus  not  surprising  that  the  single 
solitary  wave  solution  to  the  HD  is  an  implicit  one.  For  a  better  understanding  we  compare, 
wherever  possible,  different  steps  in  every  conversion  procedure.  For  instance,  it  is  realised 
that  the  Schwarzian  derivative,  which  occurs  here  in  the  context  of  Painlev4  analysis,  may 
be  conceived  of  as  a  combination  of  the  Miura  and  the  Cole-Hopf  transformations. 

We  next  turn  to  transforming  the  known  solution  (s)  of  one  equation  to  that  of  the 
other.  This  is  accomplished  in  Section  5.  Of  particular  interest  is  the  case  where,  by 
starting  from  the  cusp-soli  ton  solution  of  the  HD,  we  are  able  to  derive  a  new  closed  form, 
though  singular,  solution  of  the  KdV  equation. 


2.  Derivation  of  the  Harry  Dym  equation 

In  this  Section  we  shall  derive  the  HD  equation  by  (a)  starting  from  the  classical  string 
problem,  (b)  employing  the  Lax  operator  technique  and  (c)  tising  a  heuristic  approach. 

Method  (a): 


Consider  the  ODE 


V'l.xx  = 


r*(x;f) 


(1) 


which  models  the  classical  string  problem  where  the  string,  for  instance,  has  a  varying  elas¬ 
tic  constant  (Sabatier  1979a,  b,  1980,  Dijkhuis  and  Drohm  1979,  Yi-Shen  1982,  Calogero 
and  Degasperis  1982).  In  (1),  A  denotes  the  eigenvalue  and  r(z;t)  is  a  bounded,  positive 
function  of  z.  Furthermore,  we  assume  that  r(z;t)  — »  0  as  |z|  oo.  For  r(z;t)  =  1,  Eq. 
(1)  reduces  to  the  standard  Schrodinger  equation.  If  r(z;t)  ^  1,  Eq.  (l)  can  be  trans¬ 
formed  back  to  that  standard  equation  by  a  suitable  change  of  variables  (see  Section  4). 
Both  r  and  depend  on  the  parameter  t.  In  order  to  find  an  integral  representation  for 
the  solution  of  (1)  that  also  incorporates  the  conditions,  we  first  reduce  (1)  to  an  equivalent 
system  of  first-order  equations.  Defining  =  ^i,x>  (I)  may  then  be  rewritten  as 


Imsim 


We  remark  here  that  the  above  decomposition  is  not  unique  in  the  sense  that  the  state 
variables  and  ^*2  could  be  chosen  differently. 

It  was  noted  by  Schlesinger  (1912)  that  solvable  nonlinear  PDEs  arise  as  the  integra- 
bility  conditions  when  a  linear  ODE  as  in  (2)  is  deformed  to,  for  instance, 


».  =  Ar«,  fi) 


(3) 


in  such  a  way  as  to  preserve  some  eharaeteristie  (for  instance,  the  spectrum  A)  of  the 
equation.  Thus,  setting  At  =  0  (isospectral  case),  the  integrability  conditions 


(4) 


yield  the  structure  equation 


Mt-N^  +  [M,N]  =  0, 


(5) 


where  we  define  the  commutator  [M,  N]  =  MN  —  NM.  Incorporating  the  explicit  forms 
for  M  and  N  (from  (2b)  and  (3))  leads  to  four  coupled  equations,  one  of  which  describes 
the  evolution  of  r(z,t)  : 

c  +  4b  =  a*, 

••A 


£>  -  A  =  B*, 


-  b-c  =  d. 


(ба) 

(бб) 
(6c) 

(6d) 


From  (6a)  and  (6d),  it  follows  that 


D  =  -A, 


(7) 


where  we  have  neglected  a  constant  of  integration.  Using  (6a,b)  and  (7)  in  (6c),  the 
evolution  of  r  may  be  rewritten  as 


li  _ _ Lr  J-IiR—"®* 

"3x  "I"  " 


4A 


2  • 


(8) 


The  choice  B  =  — 4Ar  (Yi-Shen  1982)  leads  straightforwardly  to  the  HD  equation 


n  =  r-rs*. 


(9) 


Method  (b): 


Following  Lax  (1968),  we  can  derive  an  eigenvalue  problem  of  the  tyi>e  =  Xip  with 
At  =  0  for  a  nonlinear  evolution  equation  expressible  in  the  form 


Lt  —  [B,L], 

4 


(10) 


where  L  and  B  are  linear  spatial  operators  (respectively  symmetric  self-adjoint,  and  anti¬ 
symmetric)  which  may  depend  on  r  and  its  spatial  derivatives.  Now  the  ODE  (1)  may  be 
modified  to 

+  =  (11) 

where  ^  is  a  new  wave  function  defined  as  ^  and  where  m  is  an  integer.  From 

the  construnt  of  self-adjointness,  it  readily  follows  that  m  =  1,  yielding 

L  =  rgjjr,  (12) 

SO  that  (10)  becomes 

+  =  (13) 

To  construct  the  antisymmetric  operator  £,  we  set,  following  Sabatier  (1979b), 

B  =  AL  +  LA  (14) 

where  A  is  antisymmetric;  thus  [B^L]  =  [A,L]L+L[A,L].  Comparing  with  (13),  it  follows 
that 

K£1  =  7.  (15) 

provided  the  LHS  of  the  above  equation  reduces  to  a  multiplicative  operator  and  (15) 
yields  the  required  nonlinear  PDE.  Noting  that 

^  =  =  O') 

we  set  i4  =  5  -1-  T,  so  that 

(>i,L]  =  5(5,r)-(5,r]r 

=  (5.[5,r]]-i-[5,r)(5-r) 

=  »'("2x),  -  "xJ-j* 

=  r’rsx,  (17) 

where  we  have  used  the  observations  T  —  S  =  r  and  [5,r]  =  rr2i.  From  (15)  and  (17), 
the  HD  equation  (9)  follows. 

Method  (e): 

The  heuristic  approach  to  constructing  nonlinear  evolution  and  wave  equations  in¬ 
volves  deriving  the  linear  PDE  from  a  (known)  dispersion  relation  and,  thereafter,  in¬ 
corporating  nonlinear  corrections  to  the  (linear)  phase  velocity  cq  (Korpel  and  Banerjee 
1984).  Consider  the  dispersion  relation 


where  u  denotes  the  angular  frequency  and  k  the  propagation  constant.  The  linear  PDE  is 
obtained  through  replacing  u  and  k  by  their  respective  operators  u>  — »  —idfdi,  k  -*  idjdx 
in  (18)  : 

fi-\-CQfi  +  'ifzi  =  0.  (19) 

The  nonlinear  extension  is  now  done  by  assuming  that  'y,  instead  of  Cq,  is  a  nonlinear 
function  of  f 

'y  =  'y3f®,  (20) 

where  73  is  constant.  Incorporating  this  in  (19)  together  with  the  scalings 

i  =  x-Cot,  t  =  t,  r(i,t)  = -(73)*r(xJ),  (21) 

we  readily  retrieve  the  HD  equation  as  in  (9). 


3.  Derivation  of  the  implicit  single  solitary  wave  solution 

In  this  section,  we  employ  direct  integration  to  construct  the  (implicit)  one  solitary 
wave  solution  of  the  HD  equation  previously  found  by  Wadati  tt  al  (1980)  using  the  1ST. 
We  are  guided  by  the  fact  that  the  solitary  wave  solution  has  a  singularity  in  its  derivative 
at  its  (bounded)  peak  value  owing  to  the  presence  of  a  transcendental  phase  c(x,t)  (see 
figure  1).  We  therefore  seek  for  a  solution  to  (9)  in  the  form 

r(x,  <)=/’(/)  (22a) 

where 

/(x,t)  =  if(x  -  xo  +  Vt)  + /iCc(x,t)  (226) 

with 

/fc(x,t)  =  G{f).  (22c) 

In  the  above  scheme,  xq  is  a  constant,  K  is  the  wavenumber  and  V  is  the  anticipated 
velocity  pf  the  solitary  wave. 

At  this  point,  we  wish  to  speculate  on  the  physical  basis  of  implicit  solutions  of 
nonlinear  PDEs  in  general.  Readers  are  reminded  here  of  the  real  exponential  approach  to 
solving  nonlinear  evolution  and  wave  equations  (Korpel  1978,  Hereman  tt  al  1985,  1986) 
where  the  final  solution  is  assumed  to  be  built  up  from  the  nonlinear  mixings  of  real 
exponential  solutions  to  the  lintar  dispersive  part  of  the  PDE.  Alternatively,  we  may  think 
of  constructing  a  particular  solution  from  the  solution  to  the  nonlinear  nor^dispersivt  part 
of  the  PDE.  This  is  a  valid  conjecture,  since  the  nonlinear  nondispersive  part  of,  say,  the 
KdV  equation  in  tt(X,r),  i  f*  ,  the  nonlinear  kinematic  equation  (Whitham  1974), 


Ur  +  outtx  =  0, 


(23) 


possesses  shock  wave  solutions  that  are  intrinsically  implicit  : 


u{X,T)  =  g{X-au(X,T)T). 

(24) 

In  a  more  general  sense,  we  may  therefore  think  of  solutions  to  the  entire  nonlinear  dis¬ 
persive  PDE  to  be  of  the  form  (Banerjee  and  Hereman  1988a) 

t.(x,r)  =  f(/) 

(25a) 

with 

nX,T)=H,{f)X-H:,[f)T  +  H^[f). 
Examining  (22),  we  note  that  it  fits  the  scheme  described  above. 

(256) 

Returning  to  the  construction  of  the  implicit  solution  to  (9), 

that 

d  KV  d  d  K  d 

dx  ^\-G/df' 

we  first  note  from  (22) 

(26) 

Eq.  (9)  then  transforms  to 

K^F^  dP^l-GfUri-Gf^' 

(27) 

A  first  integration,  followed  by  multiplication  by  Ff  and  a  subsequent  integration  yields 


=  e,(l  -  G,fF  +  c,{l  -  G,f  +  (28) 

where  ci,C2  are  integration  constants. 


To  solve  for  F,  we  must  first  eliminate  G  from  (28)  by  assuming  a  relation  between 
Gf  and  F.  The  obvious  choice, 

l-Gf  =  F,  (29) 

reduces  (28)  to  an  ODE  which  can  be  readily  solved  either  by  employing  the  real  expo¬ 
nential  approach  (Hereman  et  al  1985,  1986)  or  by  direct  integration.  Using  the  latter,  we 
get 


2 


f  =  j  (^F+«,F’  +  c,F’)  'dF, 
where  the  integration  constant  can  be  absorbed  in  xq.  For  the  choice 


(30) 


(31) 


Eq.  (30)  may  be  readily  evaluated  as 


From  (29)  we  then  obtain 


G(/)=^Unh(||/). 


For  later  use,  we  write  down  the  implicit  solution  to  the  HD  equation  in  the  original 
variables  x  and  t  : 


r(x,t)  =  tanh'f-— (z  -  zq  +  K/  +  t(z,t))] 

£i 


c(a,0  =  ^tanh(^(z-zo  +  Vt  +  £(z,f))j. 


This  type  of  cusp-soliton  solution  (which  is  plotted  in  figure  1)  has  also  been  obtained 
for  coupled  systems  of  evolutions  equations.  Kawamoto  (1984a, b,  1985a)  discusses  two 
examples  :  an  Ito-type  system  and  a  normalised  Boussinesq  equation. 


4.  Connections  between  the  HD  and  the  KdV  equations 


Having  derived  the  cixsp-solitary  wave  solution  to  the  HD,  we  would  like  to  investigate 
how  this  solution  maps  to  particular  (hopefully  new)  solutions  of  the  KdV.  Toward  this 
goal,  we  first  reconstruct  the  transformations  between  the  HD  and  the  KdV  equations 
according  to  the  following  scheme  : 


(1)  by  deriving  the  HD  from  the  KdV  equation  using 

(i)  the  Backlxmd  transformation  obtained  by  Weiss  (1983,  1986), 

(ii)  the  Sehwarzian  transformation  proposed  by  Ibragimov  (1981,  1985); 

(2)  by  retracting  the  KdV  from  the  HD  using  the  Cole-Hopf  and  the  Miura  transforma¬ 
tions  as  suggested  by  Kawamoto  (1985b); 

(3)  by  transforming  the  spectral  problem  associated  with  the  HD  equation  to  the  one  for 
the  KdV  equation  (Calogero  and  Degasperis  1982). 


We  remark  that  in  addition  to  every  transformation  above,  we  need  an  explicit  (im¬ 
plicit)  to  implicit  (explicit)  transformation  in  going  from  the  KdV  (HD)  to  the  HD  (KdV). 
This  explains  why  the  solitary  wave  solution  of  the  HD  equation  is  implicit  in  nature. 
Detailed  calculations  on  the  conversion  of  the  cusp-soliton  solution  of  the  HD  to  particular 
solutions  of  KdV  (and  vice  versa)  will  be  presented  in  Section  5. 


4.1.  Transformation  from  the  KdV  to  HD 


(i)  Backlund  transformation  method. 

This  conversion  process  from  the  KdV  to  the  HD  can  be  summarized  in  the  following 
three  steps  : 


1.  We  first  use  the  auto-Backlund  transformation 


tt{x,T)  =  -|in*(;f.r)i,x  +  uj(x,r) 
a 


>>  • . . 


mum 


^  O  V  V  'tf  X."  V7  W"  ^  O  O 


W'/  ▼  .K  •^W  "  .  »  « 


derived  by  Weiss  (1983,  1986)  in  the  context  of  the  weak  Painlevd  analysis  of  the  KdV 
equation 

Ur  +  auux  +  usx  —  0-  (35) 

The  function  ^(X,  T)  plays  a  crucial  role  as  a  singular  or  poU  manifold  (Newell  et  al 
1987)  in  the  Painlev^  formalism  and  in  Hirota's  bilinear  method  as  a  new  dependent 
variable  (Matsuno  1984,  Gibbon  et  al  1985).  We  remark  that  U2[X^T),  like  u{X,T), 
must  also  be  a  solution  to  (35).  Hence,  from  subtracting  the  respective  equations,  we 
obtain 

(36) 


(u  —  tt2)t  +  (u  —  a2)3X  +  "2  1(“  ~  “2)^  +  2U2(u  —  U2)]x  =  0- 


We  now  substitute  for  u  —  U2  from  (34)  into  (36)  and  integrate  once  w.r.t.  x.  After  a 
little  algebra  and  equating  the  coefficients  of  4>~^  and  4>~^  we  obtain  the  two  nontrivial 
relations  : 

^xr  +  ot^2XU2  +  ^4X  =  0.  (37a) 

4>x4>t  +  «^x^«2  +  ^4>x4>zx  —  3^2X^  =  0.  (376) 

The  coefficients  of  and  identically  vanish.  Upon  elimination  of  U2  from  (37) 
and  one  integration  w.r.t.  X,  we  get 


•t>x 


+  {^;X}  =  M, 


(38) 


where  the  Schwaizian  derivatWe  (Hille  1976)  is  defined  by 


<PX  2  <Px 


(39) 


In  (38)  n  may  depend  on  T  but  is  constant  in  X.  We  recall  the  property  of  Galilean 
invariance  of  (38)  :  If  ^{X,T]n)  denotes  the  solution  to  (38),  then 


^(X,r;0)=^(X-Mr,r;M). 


(40) 


2.  This  step  involves  the  explicit  to  implicit  transformation  by  interchanging  the  depen¬ 
dent  variable  ^  and  the  independent  variable  X.  This  is  effected  by  first  recalling  a 
property  of  the  Schwarzian  derivative  (Hille  1976),  viz.. 


(41) 


Together  with  the  observations 


^x  =  x;-*  T  constant. 


(42) 


and 


=  ^xXt  = 


Xt 


X^' 


(43) 
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m 

i 

3 


.3 


r, 


1 

\  (38)  may  be  recast  into  the  form 

- -  j 

1 

1  Xr  =  /i+^{X;<^}, 

(44) 

) 

1 

S  with  A’  =  X(^,r). 

1 

3.  Finally,  we  define  a  new  dependent  variable 

1 

J 

* 

i 

(45) 

^  and  study  the  evolution  of  w{(f>^T).  We  note  that 

- 

1  X^T  2  V 

■  wt  =  2  ~  ^ 

(46)  1 

y  which  can  be  reexpressed,  using  (44),  as 

1  WT  =  -  to^(u;u;2^  + 

S  =  w^vjz<f,  - 

(47)  : 

Consider  the  special  case  where  =  0.  Upon  defining  new  variables 

: 

. 

■ 

II 

H 

u 

II 

II 

H 

(48) 

V 

% 

^  we  readily  obtain  the  HD  equation  (9). 

N 

V 

>1 

v* 

(ii)  The  Schwarzian  transformation  method. 

"1 

||  The  next  conversion  method  from  the  KdV  to  the  HD  is  partly  motivated  by  the 

inverse  transformation  due  to  Kawamoto  (1985b).  In  essence,  Kawamoto’s  procedure 
^  involves  an  implicit  to  explicit  transformation  to  reduce  the  HD  as  in  (9)  to  an  in- 

!*¥  termediate  equation  similar  to  (38),  and  then  using  the  Cole-Hopf  transformation  on 

4>x  to  yield  the  mKdV  equation.  The  Miura  transformation  finally  links  the  mKdV 

»  to  the  KdV. 

t;-: 

To  go  from  the  KdV  to  the  HD  we  will  apply  the  above  transformations  in  reverse  J 

order  to  yield  the  intermediate  equation  (38).  In  that  sense,  this  method  differs  from  ^ 

(i)  in  the  first  step.  We  first  use  the  Miura  transfornoation  to  v{X,  T)  :  * 

1  ¥  ■  /6 
t;  “  =  ”  +‘V5'’^' 

(49a)  : 

followed  by  the  Cole-Hopf  transformation  on  the  potential  <l>x{X,T)  : 

n?  .  <f>2x 

1  S  2a  4>x' 

(496)  ; 

J  10 

i 

2 

■ 

1 

•jVh  "fc  ' 

After  some  tedious  algebra,  which  has  been  performed  using  MACSYMA,  we  can 
identify  the  relation 

f'l  (^T  +  =  0  (50) 

where  the  operator  Li  is  given  by 


I  -  ^  ^  f4>zx  o^ax^v  d  . 

'  a^4>xBX^  <t>x'^dX^  <t>x^^dX^' 


(51) 


from  which  we  recover  the  intermediate  equation  (38).  With  MACSYMA,  we  have 
obtained  yet  another  but  equivalent  representation  of  (50)  : 


ij(|^  +  W:  A-})  =  0, 


(52) 


with 

The  latter  operator  is  similar  to  the  anti-symmetric  Lax  operator  for  the  KdV,  wherein 
we  have  formally  replaced  u{X,T)  by  a  multiple  of 


Note  that  the  successive  transformations  defined  in  (49)  may  indeed  be  combined  to 
give 

«  =  |{^:A},  (54) 

which  is  the  Sehwarzian  transformation  proposed  by  Ibragimov  (1981,  1985).  The 
argument  presented  in  Ibragimov’s  book  may  be  summarised  as  follows  :  Suppose 
there  exists  a  transformation  from  (38)  to  the  form 

’*'r  +  ^3X  +  ^  (’i',  '*'x)  =  0  (55) 

having  a  nontrivial  Lie-Backlund  algebra.  Then  (38)  is  equivalent  to  either  a  linear 
equation  with  constant  coefficients  or  the  KdV  equation.  In  our  case,  (49b)  provides 
the  relevant  transformation  to  the  mKdV  equation  where  'i'  =  v;  $  =  pv^vx,  which 
is  equivalent  to  the  KdV  through  the  Miura  transformation  (49a). 


The  subsequent  steps  for  the  transformation  are  identical  to  steps  2  and  3  in  (i). 


4.S.  Transformation  from  the  HD  to  KdV 


Along  the  lines  of  Kawamoto  (1985b),  we  convert  the  HD  equation  to  the  KdV  through 
the  following  three  steps  : 

1.  We  first  employ  an  implicit  to  explicit  transformation  by  defining 


X  = 


-f, 


(56) 
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with  R[XyT)  =  r{x[X,T)^t{X,T))  representing  the  new  trainsformed  dependent  vari¬ 
able.  We  observe  that 


d  5  /  1  2\  ^ 

S  =  -ar  + 

_a_ 

dT  ^  iE2 

where  use  has  been  made  of  the  HD  equation  (9)  to  replace  all  time  derivatives  of  r 
in  terms  of  spatial  derivatives,  and  that 

d  Id 
Tt  =  Rax- 

After  a  little  algebra,  Eq.  (9)  can  be  expressed  as 


Rt  + 


R3xR^-ZR2xRxR+^Rl 


^  =0. 


In  the  above  derivation  we  have  used  the  fact  that  r(x,  t)  and  its  spatial  derivatives 
tend  to  zero  as  |x|  —*  oo. 

Before  proceeding  any  further,  we  make  the  following  remarks  : 

Remark  1  :  The  substitutions  in  (56)  are  equivalent  to 


r**  ds 

X=  /  -7— r,  T  =  t. 

Jz 


Remark  2 :  Comparing  (58)  with  (38),  it  may  be  readily  recognized  that 

R  =  (f>X’ 


2.  We  now  write 


and  find  vr(-X^»?’)  with  the  use  of  (58).  It  is  easy  to  check  that  vt  can  be  expressed 
entirely  in  terms  of  v  and  its  spatial  derivatives  through  the  mKdV  equation 

3 

wt  -  +  vzx  =  0.  (62) 


3.  The  Miura  transformation 


now  relates  (62)  to  the  KdV  equation  (35).  This  may  be  readily  verified  by  computing 
Ur  from  (63),  using  (62)  to  replace  the  time  derivative(s)  of  v  in  terms  of  v  and  its 
spatial  derivatives,  and  finally  reexpressing  in  terms  of  u. 


4.S.  Connection  between  the  respective  eigenvalue  problems 


We  will  complete  the  connection  between  the  HD  and  the  KdV  equations  by  pointing 
out  the  link  between  their  respective  eigenvalue  problems.  To  this  end,  we  first  employ 
the  implicit  to  explicit  transformation  described  by  (56),  which  is  similar  to  the  one 
employed  by  Calogero  and  Degasperis  (1982).  The  eigenvalue  problem  (1)  for  the  HD 
equation  is  then  equivalent  to 


-  ^xRx  _  X  a 

- R - 


where  0[X]T)  =  V'i(a:(X,r);t(r)).  Now  defining 

(1{X;T)  =  R-h{X;T) 


and  eliminating  6  from  (64),  we  obtain 


(66a) 


or,  using  (60), 


^2X  + 


Note  that  if  we  adopt  the  Cole-Hopf  transformation  (61),  then  (66a)  may  be  also 


written  as 


n2X  +  (7(v*-2t;A:)-A]n  =  0 

4 


and  thence  as 


^2x  +  (®  ~  A]n  —  0, 

provided  u(A',T)  is  a  solution  to  the  KdV  equation 


Ur  —  6ttttx  +  ttsx  =  0. 


Thus,  we  have  transformed  the  eigenvalue  problem  for  the  HD  to  that  for  the  KdV, 
where  the  corresponding  eigenfunctions  ^i(=  9)  and  fi  are  related  through  (65).  The 
explicit  form  for  R{X,T)  will  be  discussed  in  the  following  section. 


5.  Connections  between  the  solutions  of  the  HD  and  KdV  equations 


Using  the  knowledge  of  the  connections  between  the  HD  and  the  KdV  equations 
discussed  in  the  previous  section,  we  shall  proceed  to  transform  the  known  solution  (s)  of 
each  equation  to  the  solution  (s)  of  the  other.  The  transformation  of  the  implicit  solution 
of  the  HD  to  explicit  solution(s)  of  the  KdV  is  more  interesting  and  will  be  taken  up  first. 


We  shall  first  apply  the  implicit  to  explicit  transformation  prescribed  by  (56)  to  the 
cusp-solitary  wave  solution  of  the  HD  as  in  (33).  From  the  latter  equation 


hence,  from  (56), 


-  =  1  +  «x  , 

r 


X  =  x-t-e(x,t) 


since  x  +  Vt  +  e(x,t),  which  is  the  argument  of  the  tanh^  function  in  (33a),  tends  to  zero 
as  \x-\-Vt\  -*  oo  (see  figiire  1).  Consequently,  from  (33a)  and  the  definition  of  /2  as  in 
(56), 

R{X,T)  =  tanh"[^(A'-  Xo  -  VT)\,  (72) 

where  Xo(=  Xq)  is  a  constant. 

For  later  use,  we  will  write  down  the  expression  for  ^  from  (60).  Using  (72),  and  one 
integration  w.r.t.  X,  we  get 

^  =  X-X,(T)-  ^tuthl^(X  -Xo-  vr)l,  (73) 

where  Xi(T)  is  an  integration  constant. 

From  (61)  and  (72),  we  arrive  at  the  singular  solution  to  the  mKdV  equation  (62)  : 

v(X,  r)  =  2\/Kco8ech(\/V(X  -  Xo  -  VT)].  (74) 

The  solution  to  the  KdV  equation  (35)  may  now  be  calculated  using  the  Miura  transfor¬ 
mation  (63).  Straightforward  algebra  yields  the  familiar  soliton  solution 

ti{X,T)  =  — -  Jfo  -  vr)).  (75) 

a  z 

An  alternate  solution  of  the  KdV  may  be  derived  by  starting  from  (73)  and  employing 
the  auto-Backlund  representation  of  u  in  terms  of  as  in  (34).  Indeed,  since  <f>  should  also 
satisfy  (38),  it  is  readily  checked  upon  substituting  (73),  that  Xi{T)  =  constant  =  Xi  and 
fi  =  0.  Now,  U2  can  be  computed  using  either  of  the  relations  (37a,b).  This  gives  again 
the  one-solitary  wave  solution  for  uj,  as  in  (75).  Returning  to  (34),  a  new,  but  singular, 
solution  to  the  KdV  equation  reads 

’  “  (tanh|i^(X-Yo-Vr)|-^(Y-Y,))7  '  ' 


We  remark  that  the  inverse  process  of  transforming  the  familiar  sech^  solution  of 
the  KdV  to  that  of  the  HD  leads  to,  for  instance,  a  rational  solution  r(x,t)  =  ^(x  —  l). 


Other  possible  solutions  to  the  HD  and  the  KdV  equations  are  presently  under  investigation 
(Banerjee  and  Hereman  1988b). 

6.  Conclusions 

In  retrospect,  we  have  provided  some  straightforward  derivations  of  the  HD  equation. 
A  direct  integration  method  was  used  to  derive  a  simple  particular  implicit  solution.  We 
have  consolidated  the  various  links  between  the  HD,  the  KdV  and  the  mKdV  equations, 
and  used  these  to  provide  connections  between  their  solutions. 

Although  at  this  stage,  the  HD  equation  is  more  of  theoretical  significance  than  of  ap¬ 
plicative  relevance,  we  hope  that  this  paper  contributes  not  only  to  provide  a  better  insight 
into  cusp  solitary  waves,  but  also  serves  as  a  review  on  the  HD  equation.  The  straight¬ 
forward  mathematical  techniques  employed  throughout  this  paper  should  prove  helpful  in 
deriving  implicit  solutions  to  other  nonlinear  PDEs  and  in  (as  a  spin-off)  discovering  new 
solutions  to  the  KdV  equation.  Work  on  this  is  currently  in  progress. 
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